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Recent proton spin-lattice relaxation-time (Ti) measurements on the ferrimagnetic chain com- 
pound NiCu(C7H6N206)(H20)3-2H20 are explained by an elaborately modified spin-wave theory. 
We give a strong evidence of the major contribution to l/Ti being made by the three-magnon scat- 
tering rather than the Raman one. 



PACS numbers; 75.10.Jm, 75.50.Gg, 76.50.+g 

Low-frequency spin dynamics in magnetic systems is 
a long-standing problem and has recently attracted re- 
newed interest due to the significant progress in design- 
ing low-dimensional materials such as chains and ladders. 
Nuclear magnetic resonance (NMR) is a powerful probe 
to their dynamic properties. The nuclear spin-lattice re- 
laxation time Ti is in particular eloquent of the collective 
motions of electronic spins and therefore we take a great 
interest in microscopically interpreting it. The spin- wave 
formalism has played a crucial role in this context. Van 
Kranendonk, Bloom [1], and Moriya [2] made their pio- 
neering attempts to describe 1/Ti in terms of spin waves. 
Oguchi and Keffer [3] further developed the spin-wave 
analysis considering the three-magnon nuclear relaxation 
mechanism as well as the Raman one, whereas Pincus and 
Beeman [4] claimed that the three-magnon process was 
considerably underestimated in their argument, revealing 
further relaxation mechanism. 

The spin-wave excitation energy is usually much larger 
than the nuclear resonance frequency and thus the single- 
magnon relaxation process is rarely of significance. The 
Raman process consequently plays a leading role in the 
nuclear spin-lattice relaxation. Because of the (45')^^- 
damping factor to the Holstein-Primakoff magnon series 
expansion, the multi-magnon scattering is much less con- 
tributive within the first-order mechanism, where a nu- 
clear spin directly interacts with spin waves through the 
hyperfine coupling. However, the second-order mecha- 
nism, where a nuclear spin flip induces virtual spin waves 
which are then scattered thermally via the four-magnon 
exchange interaction, may generally enhance the relax- 
ation rate. The Pincus-Beeman spin-wave nuclear relax- 
ation theory is thus fascinating but works only far below 
the three-dimensional transition temperature. The con- 
ventional spin-wave theory applied to low-dimensional 
magnets ends in failure with diverging magnetizations. 
In such circumstances, Takahashi [5] gave a fine descrip- 
tion of the low-dimensional ferromagnetic thermodynam- 
ics at low temperatures in terms of modified spin waves. 
His idea of introducing a constraint on the magnetiza- 
tion was developed for antiferromagnets [6,7] and fer- 
rimagnets [8,9]. Random-bond ferromagnets [10] and 
frustrated antiferromagnets [11-13] were also discussed 
within this renewed spin-wave scheme. 



The ferrimagnetic modified spin- wave theory is partic- 
ularly powerful to investigate both static [14,15] and dy- 
namic [16,17] properties. One-dimensional ferrimagnets 
have lately attracted much attention especially in the 
context of designing molecule-based ferromagnets [18]. 
A series of bimetallic chain compounds [19] are typical 
examples and some of them were indeed assembled into 
a ferromagnetic lattice. Another approach [20] consists 
of bringing into interaction metal ions and stable organic 
radicals. Homometallic systems can exhibit distinct fer- 
rimagnetism of topological origin [21,22]. Such syn- 
thetic endeavors have stimulated several experimentalists 
[23,24] to measure Ti on ferrimagnetic chain compounds. 
Thus motivated, here we make a systematic spin-wave 
analysis on the nuclear spin dynamics in one-dimensional 
Heisenberg ferrimagnets, which has been pending for the 
past decades without any suitable spin-wave scheme. Our 
goal is to show a strong evidence of the proton spin relax- 
ation in the title compound being mediated by the three- 
magnon scattering rather than the Raman one. 

First of all our scheme [8] of modifying the spin-wave 
theory is distinct from the original idea proposed by 
Takahashi [6] and Hirsch et al. [7]. Their way of sup- 
pressing the divergent sublattice magnetizations consists 
of diagonalizing an effective Hamiltonian with a Lagrange 
multiplier included subject to zero staggered magnetiza- 
tion. The thus-obtained energy spectrum necessarily de- 
pends on temperature and fails to reproduce the Schot- 
tky peak of the specific heat [25]. In order to obtain 
better thermodynamics, we first diagonalize the Hamil- 
tonian keeping the dispersion relations free from temper- 
ature and then introduce a Lagrange multiplier in order 
to minimize the free energy subject to zero staggered 
magnetization. This scheme is highly successful in de- 
scribing the magnetic susceptibility as well as the specific 
heat [25] and therefore guarantees our exploration of the 
one-dimensional ferrimagnetic dynamics over a wide tem- 
perature range. 

We consider ferrimagnetic Heisenberg chains of alter- 
nating spins S and s, as described by the Hamiltonian 

TV 

n=Y, [jSn ■ {Sn-1 + Sn) ~ {gsS'^ + 5s<jMBi?] ■ (1) 
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Introducing bosonic operators for the spin deviation in 
each sublattice via — (25— aL:ja-:i)^^^a_:,;, — S— 

and assuming that 0[S) = 0{s), we expand the Hamilto- 
nian with respect to 1/5 as = H2+'Hi+Ho + 0{S~^), 
where Hi contains the 0{S'^) terms. H2 = —2SsJN is 
the classical ground-state energy, while Hi describes lin- 
ear spin-wave excitations and is diagonalized in the mo- 
mentum space as 

Hi^~iS + s)JN- [gs{S+^)-g,is + ^)]^iBHN 
S + s + {gs-gs)fiBH/2J 
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where a^.j. = a^.j,cosh0fe + a-cs-ksmhOk^ provided 
tanh26'fe = 2V5s cos(fc/2)/[5 + s + {gs - gs)nBH/2J], 
creates a spin wave of ferromagnetic {a — — ) or an- 
tiferromagnetic (a = +) aspect [26], whose energy is 
given by cj^ — uj^ + (j[S — s — {gs + gs)fJ,BH/2J] with 
iOk = {[S + s + {gs - g,)fiBH/2J]^ - ASscos^k/2)y/\ 



Minimizing the free energy under the condition of zero 
staggered magnetization [16], we obtain the optimum dis- 



tribution functions 
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The hyperfine interaction is generally expressed as 

+5s/iB?i7N/+E„(5^« + KO, (3) 



where A"^ {B^) is the dipolar coupling tensor between the 
nuclear and nth larger (smaller) electronic spins. Since 
Ho and Hht are both much smaller than Hi, they act as 
perturbative interactions to the linear spin-wave system. 
If we consider up to the second-order perturbation with 
respect to V = Hq + Hm, the probability of a nuclear 
spin being scattered from the state oi — m to that of 
— m + \ is given by 
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(a) Direct process 
(first order) 



(b) Raman process 
(first order) 
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(c) 3-magnon process 

(second order) 
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FIG. 1. Diagrammatic representation of various nuclear spin-lattice relaxation processes. Solid arrows, designating spin 
waves which are emitted in the first-order processes, induce a nuclear spin flip (x) via the hyperfine interaction, while broken 
arrows, depicting four-magnon exchange correlations, thermally scatter the first-order spin waves as virtual excitations, where 
spin waves of ferromagnetic and antiferromagnetic aspect are distinguishably drawn by straight and wavy arrows, respectively, 
(a) The first-order direct (single-magnon) relaxation processes; (b) The first-order Raman (two-magnon) relaxation processes; 
(c) The first-order and second-order three-magnon relaxation processes, where q = —k^ = ki — k2 — k^, are related to each other 
through nonhnear equations and are therefore inseparable. Considering the nuclear-electronic energy conservation, processes in 
solid and dotted frames are of great and little significance, respectively, whereas those in broken frames are relevant according 
to the constituent spins S and s. 
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where i and / designate the initial and final states of 
the unperturbed electronic-nuclear spin system. Then 
we find that Ti = (/ - m)(/ + m + 1)/2W. Equation 
(4) contains various relaxation processes but their ex- 
plicit formulae will be presented elsewhere. We instead 
diagrammatically show them in Fig. 1. Due to the con- 
siderable difference between the nuclear and electronic 
energy scales, htujq <C J, the direct process, involving 
a single spin wave, is rarely of significance. Consider- 
ing further that the antiferromagnetic spin waves are 
higher in energy than the ferromagnetic ones, w^T < w^, 
at moderate fields, the intraband spin- wave scattering 
dominates the Raman relaxation rate l/T^ , whereas 
both the intraband and interband spin-wave scatterings 

fr>\ 

contribute to the three-magnon relaxation rate . 
Within the first-order mechanism, 1 /T^ is much smaller 
than [3]. However, the first-order relaxation rate is 

generally enhanced through the second-order mechanism. 
We consider the leading second-order process, that is, the 
exchange-scattering-induced three-magnon relaxation, as 
well as the first-order process. The second-order single- 
magnon and Raman processes, containing three and two 
virtual magnons, respectively, are much more accidental 




\ \ ksT/J = 


0.1 




Raman 




- 3-magnon 




-total 


\ k^T/J = 


= 0.5 




1.0 





10-' 10-5 10-4 

FIG. 2. Modified spin-wave calculations of typical temper- 
ature (the left three) and field (the right three) dependences 
of the nuclear spin-lattice relaxation rate, where gs ~ Qs = g, 
A^B^ = 1, and {B-/B'f = 4. and l/rf> 

are plotted by dotted and broken lines, respectively, while 
+ l/Tj^^' = 1/Ti, which is observable, by solid lines. 
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FIG. 3. The crossover point as a function of temperature 
and an applied field, where gs = gs and (_B~/_B^)^ — 4. 



due to the momentum conservation and much less 
contributive due to the (45)~^-damping factor in the 
Holstein-Primakoff magnon series expansion. 

We calculate the case of {S,s) = (1,5)7 which is rele- 
vant to several major materials [20,22,27], assuming that 
the Fourier components of the coupling constants have 
little momentum dependence [24] as J^n e''^"^^ = ~ 



and J2r. 
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shows 1 /Ti as a function of temperature and an applied 
field. The exchange-scattering-enhanced three-magnon 
relaxation rate generally grows into a major contribution 
to 1/Ti with increasing temperature and decreasing field. 
As temperature increases, decreases at fc ~ but 
otherwise increases [28]. In one dimension, excitations 
at fc ~ predominate in the Raman process, while all 
the excitations are effective in the three-magnon process. 

f2l (3) 

l/T\^> and l/T\^> are hence decreasing and increasing 
functions of temperature, respectively, unless tempera- 
ture is so high as to activate the antiferromagnetic spin 



waves. The field dependences of l/Tp' and l/T^'^'' are 
also in striking contrast. At moderately low tempera- 
tures and weak fields, fwj^ <C k^T <C J, we find that 
l/rf ^ cx e~^as+g.)iJ.i,H/2k^T j^^ij^^^^ i2}.^T), where Kq is 
the modified Bessel function of the second kind and be- 
haves as ii:o(?iwN/2fcBT) ~ 0.80908-ln(nwN/fcBT). Thus 

(2) 

the field dependence of 1 /Tf is initially logarithmic and 
then turns exponential with increasing field. Equation 
(4) claims that 1/T{ is much less analyzable but should 
exhibit much stronger power-law diverging behavior with 
decreasing field. Therefore, the three-magnon relaxation 
process predominates over the Raman one at weak fields. 
In Fig. 3 we plot the crossover points on which 

f21 f3) 

1/T{ = l/Tf . A Raman-to-three-magnon crossover 
may generally be detected with increasing temperature 
and decreasing field. The ferrimagnetic nuclear spin- 
latticerelaxation is sensitive to another adjustable pa- 
rameter /B'^, that is, the location of the probe nuclei. 
At the special location of A''/B^ ~ (ds/ds f ^ {S/s)", 
where ds (ds) is the distance between the nuclear and 
larger (smaller) electronic spins, the a excitation mode 
hardly mediates the nuclear spin relaxation [28]. For 
{S,s) = (1, i), the lower-lying ferromagnetic spin waves 
are almost invisible to the nuclear spin located as /B'^ 
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FIG. 4. Froton spin relaxation-time measurements on 
NiCu(C7H6N206)(H20)3-2H20 [24] compared with our the- 
ory, (a) 1/Ti as a function of temperature at various values 
of an apphed field; (b) 1/Ti as a function of an apphed field 
at 280 K, where I/t}^' and are also plotted by dotted 

and broken lines, respectively. 
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~ 1/2 and therefore its relaxation rate stays extremely 
small. Any Ti measurements should be performed away 
from such magic points. 

We are further excited to compare our theory with re- 
cent experimental findings. Fujiwara and Hagiwara [24] 
measured Ti for proton nuclei in the bimetallic chain 
compound NiCu(C7H6N206)(H20)3-2H20 [27] compris- 
ing ferrimagnetic chains with alternating octahedral Ni^"*" 
and square- pyramidal Cu^"*" ions. The measured sus- 
ceptibility [29] suggests that J/kn ^ 121 K, gs = 2.22, 
and gs = 2.09. Comparative measurements on the D2O- 
substituted samples [24] indicate that protons relevant to 
the Ti findings are located in close vicinity to Cu spins. 
Then, considering that the dipolar coupling strength is in 
proportion to the inverse cubic distance, we may set the 
coupling constants for / B'^ = 0. Conditioning further 
that = 1.37 X IO^^tVJ and {B-fB'')'^ = 5, which 
can be consistent with the crystalline structure [27], wc 
plot calculations together with the observations in Fig. 4. 
The ferromagnetic and antiferromagnetic spin waves con- 
tribute different temperature dependences to 1/Ti and 
give the decreasing and then increasing behavior in Fig. 
4(a). Considering that there may be larger uncertainty 
in the experimental findings for 1/Ti at lower tempera- 
tures and weaker fields [24], the theoretical and exper- 
imental findings arc in good agreement and the slight 
discrepancy between them may be attributable, for in- 
stance, to weak momentum dependence of Bl and the 
protons of wide distribution. Figure 4(b) more impres- 
sively demonstrates the relevance of the three-magnon 
scattering to the proton spin relaxation. The strong field 
dependence can never be explained by the Raman process. 
Since 1 /r| within the first-order mechanism stays much 
smaller than the observations, the exchange-scattering- 
induced three-magnon process is essential in interpret- 
ing such accelerated relaxation. We are eager to have 
reliable observations at lower temperatures and weaker 
fields. We call for more extensive NMR measurements 
using as probes ^H, ^^Cu, and ^^Mn nuclei on the family 
material MnCu(C7H6N206)(H20)3-2H20 [27] as well as 
that of present interest. 

There exist pioneering Ti measurements on the lay- 
ered ferromagnet CrCla [30] and the coupled-chain an- 
tiferromagnet CsMnCl3-2H20 [31], which give evidence 
of the relevant three-magnon scattering. However, they 
are both, in some sense, classical findings under the ex- 
isting three-dimensional long-range order. Without any 
reliable spin-wave formulation in one dimension, no au- 
thor has explored quantum, ferrimagnetic dynamics with 
particular interest in multi-magnon scattering beyond 
the Raman mechanism. The present calculation is the 
first evidence of the three-magnon scattering dominating 
the one- dimensional nuclear spin relaxation and motivate 
extensive T\ measurements on various one-dimensional 
quantum ferrimagnets. 
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